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' Abstract — In a recent review paper [Phys. Reports 214 (1992) 339] we proposed, within con- 

■ ventional quantum mechanics, new definitions for the sub-barrier tunnelling and reflection times. 

^ , Aims of the present paper are: (i) presenting and analysing the results of various numerical 

^1 calculations (based on our equations) on the penetration and return times < rpen >, < TRet >) 

^ ■ during tunnelling inside a rectangular potential barrier, for various penetration depths Xf; (ii) 

^ , putting forth and discussing suitable definitions, besides of the mean values, also of the variances 

^ [ (or dispersions) D tt and D tr for the time durations of transmission and refiection processes; 

■ - - ' (iii) mentioning, moreover, that our definition < tt > for the average transmission time results 

to constitute an improvement of the ordinary dwell-time r^™ formula: (iv) commenting, at 
last, on the basis of our new numerical results, upon some recent criticism by C.R. Leavens. 
We stress that our numerical evaluations confirm that our approach implied, and implies, the 
existence of the Hartman effect: an effect that in these days (due to the theoretical connections 
between tunnelling and evanescent-^wave propagation) is receiving — at Cologne, Berkeley, Flo- 
rence and Vienna — indirect, but quite interesting, experimental verifications. Eventually, we 
briefly analyze some other definitions of tunnelling times. 
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1. — Introduction 

In our review articlel^l [Phys.Rep. 214 (1992) 339] we put forth an analysis of the main theo- 
retical definitions of the sub-barrier tunnelling and reflection times, and proposed new definitions 
for such durations which seem to be self-consistent within conventional quantum mechanics.*^ 
In particular, the "prediction" by our theory!^! of the reality of the Hartman effect^"^^ in tun- 
nelling processes has recently received — due to the analogyl^l between tunnelling electrons and 
evanescent waves — quite interesting, even if indirect, experimental verifications at Cologne, 
Berkeley, [^1 Florence'^l and Vienna J^l 

Main aims of the present paper are: presenting and analysing the results of several numeri- 
cal calculations of the penetration and return times inside a rectangular potential barrier during 
tunnelling (Sect. 3); and proposing new suitable formulae for the distribution variances of the 
transmission and reflection times (Sect. 2). The results of our numerical evaluations seem 
to confirm that our approach is physically acceptable, and that it moreover implied, and implies, 
the existence of the so-called "Hartman effect " even for non-quasi-monochromatic packets. 

The present research field, however, is developing so rapidly and in such a controversial 
manner, that it may be convenient to add here — before all — some brief comments about a few 
papers appeared during the last three years: 

(i) First, let us mention that we had overlooked a new expression for the dwell-time r^^ 
derived by Jaworsky and Wardlaw^^'^l 

TD-(xi,Xf;A:) = [j^^dtt J{x,,t) - J{x,,t)) [^^^dt J^xut))'^ , (1) 

which is indeed equivalents^] to our eq.(16) of ref.[l] (all notations being defined therein): 

r^-{x,,xf,k) = (jZ^dtj:^^ dx p{x,t)) (jZ^dt Mxut)y' . (2) 

This equivalence reduces the difference, between our definition < tt > of the average transmis- 
sion time — under our assumptions — and quantity r^™, to the difference between the average 
made by using the positive-definite probability density dt J+(x, t)/ J^^dt J+{x,t) and the av- 
erage made by using the ordinary "probability density" dtJ{x,t)/ J^^dt J{x,t). Generally 
speaking, the last expression is not always positive definite, as it was explained at page 350 of 
ref.[l], and hence does not possess any direct physical meaning. 

(ii) In ref. [9] an attempt was made to analyze the evolution of the wave packet mean posi- 
tion < x(t) > ("center of gravity"), averaged over pdx, during its tunnelling through a potential 
barrier. Let us here observe that the conclusion to be found therein, about the absence of a 

Let us take advantage of the present opportunity for pointing out that a misprint entered our 
eq.(lO) in ref.[l], whose last term ka ought to be eliminated. Moreover, due to an editorial error, in the 
footnote at page 32 of our ref. [16] the dependence of G on Ak disappeared, whilst in that paper we had 
assumed G{k -k) = C exp[-(fc - kf/{Akf]. 
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causal relation between the incident space centroid and its transmitted equivalent, holds only 
when the contribution coming from the barrier region to the space integral is negligible. 



(iii) Let us also add that in ref.[10] it was analyzed the distribution of the transmission time 
tt in a rather sophisticated way, which is very similar to the dwell time approach, however with 
an artificial, abrupt switching on of the initial wave packet. Wc arc going to propose, on the 
contrary, and in analogy with our eqs.(30)-(31) in ref.[l], the following expressions, as physically 
adequate definitions for the variances (or dispersions) Dtt and Dtr of the transmission and 
reflection time [see Sect. 2], respectively: 



and 



where 



Btt = Bt+{xf) + Bt+{xi) (3) 



Btk = Di_(xi)+Dt_(xi) , (4) 



. I-o.dtt^J±i-,t) _ lZ.dttj^ix,t) y 

Equations (3)-(5) are based on the formalism expounded in ref.[ll], as well as on our defini- 
tions for J±{x,t) in rcf.[l]. Of course, we are supposing that the integrations over J+(xf)dt, 
J-f(xi)dt and J-{xi)dt are independent of one another. We shall devote Sect. 2, below, to 
these problems, i.e., to the problem of suitably defining mean values and variances of durations, 
for various transmission and reflection processes during tunnelling. 

(iv) Below, in Sect. 4, we shall briefly re-analyse some other definitions of tunnelling dura- 
tions. 

Before going on, let us recall that several reasons "justify" the existence of different ap- 
proaches to the definition of tunnelling times: (a) the problem of defining tunnelling durations 
is closely connected with that of defining a time operator, i.e., of introducing time as a (non- 
selfadjoint) quantum mechanical observable, and subsequently of adopting a general definition 
for collision durations in quantum mechanics. Such preliminary problems did receive some clar- 
ification in recent times (see, for example, rcf.[l] and citations [8] and [22] therein); (b) the 
motion of a particle tunnelling inside a potential barrier is a purely quantum phenomenon, de- 
void of any classical, intuitive limit; (c) the various theoretical approaches may differ in the 
choice of the boundary conditions or in the modelling of the experimental situations. 
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2. — Mean values and Variances for various Penetration and Return Times during 
tunnelling 



In our previous papers, we proposed for the transmission and reflection times some formulae 
which imply — as functions of the penetration depth — integrations over time of J+{x,t) and 
J-{x,t), respectively. Let us recall that the total flux J{x,t) inside a barrier consists of two 
components, J+ and J_, associated with motion along the positive and the negative x-direction, 
respectively. Work in similar directions did recently appear in ref.[12]. 

Let us refer ourselves — here — to tunnelling and reflection processes of a particle by a 
potential barrier, confining ourselves to one space dimension. Namely, let us study the evolu- 
tion of a wave packet ^'(a;, i), starting from the initial state 5'in(a;,t); and follow the notation 
introduced in ref.[l]. In the case of uni-directional motions it is already known'^^l that the flux 
density J{x,t) = Re[{ih/m)'^{x,t) d'^*{x,t)/dx] can be actually interpreted as the proba- 
bility that the particle (wave packet) passes through position x during a unitary time-interval 
centered at t, as it easily follows from the continuity equation and from the fact that quantity 
p{x,t) = |^'(a;, i)p is the probability density for our "particle" to be located, at time t, inside a 
unitary space-interval centered at x. Thus, in order to determine the mean instant at which a 
moving wave packet ^{x,t) passes through position x, wc have to take the average of the time 
variable t with respect to the weight w{x, t) = J{x, t) / J{x, t) dt. 

However, if the motion direction can vary, then quantity w{x,t) is no longer positive defi- 
nite, and moreover is not bounded, because of the variability of the J{x, t) sign. In such a case, 
one can introduce the two weights: 



where J+{x,t) and J-{x,t) represent the positive and negative parts of J{x,t), respectively, 
which are bounded, positive definite functions, normalized to 1. Let us show that, from the 
ordinary probabilistic interpretation of p{x,t) and from the well-known continuity equation 



it follows also in this (more general) case that quantities and represented by eqs.(6), 
(7), can be regarded as the probabilities that our "particle" passes through position x during a 
unit time-interval centered at t (in the case of forward and backward motion, respectively). 

Actually, for those time intervals for which J = J+ or J = J_, one can rewrite eq.(8) as 
follows: 




(6) 




(7) 



dp{x,t) dJ{x,t) 
dt dx 



= , 



(8) 



dp>{x,t) 

at 



dJ+{x,t) 



dx 



(9.a) 
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dp<{x,t) _ dJ-{x,t) 

dt ~ dx ' ^ ^ ^ 

respectively. Relations (9. a) and (9.b) can be considered as formal definitions of dp^/dt and 
dp</dt. Let us now integrate eqs.(9.a), (9.b) over time from — oo to t\ we obtain: 

, , /■* dJ+(x,t') , , , , 

p>(x,t) = - J ' dt' (lO.a) 

, , /■* dJ-(x,t') , , , 
p< (x, t) = - dt' (10.6) 

with the initial conditions /9>(x, —oo) = p<(x, —oo) = 0. Then, let us introduce the quantities 
N>{x,oo;t) = / p>{x',t)dx'= / J+{x,t')dt' >0 (ll.a) 

Jx J —oo 

N<{-oo,x;t)= [ p<{x' ,t)dx' = - f J-{x,t')dt' >0, (11.6) 

J —oo J —oo 

which have the meaning of probabilities for our "particle" to be located at time t on the semi-axis 
(a;, oo) or (— oo, x) respectively, as functions of the flux densities J+{x,t) or J-{x,t), provided 
that the normalization condition p{x,t)dx = 1 is fulfilled. The r.h.s.'s of eqs.(ll.a) and 

(11. b) have been obtained by integrating the r.h.s.'s of eqs.(lO.a) and (10. b) and adopting the 
boundary conditions J+(— oo,i) = J_(— oo,t) = 0. Now, by differenciating eqs.(ll.a) and 
(ll.b) with respect to t, one obtains: 

dN^(x,oo,t) / ^ ^ X 
' =J+{x,t)>Q (12.a) 



^^^^%p^ = -J-(-,t)>0. (12.6) 

Finally, from eqs.(ll.a), (ll.b), (12. a) and (12.b), one can infer that: 

dNy{x,Qo\t)/dt 

w+{x,t) = ——^ (13.a) 

^ ^ A^>(a;, -oo,oo) ^ ' 

dN<{x,-<x;t)/dt 

w-(x,t) = —rr-, r— , (13.6) 

^ ' A^<(-oo,x;oo) ^ ' 

which justify the abovementioned probabilistic interpretation of w-\-{x,t) and w-{x,t). Let us 
notice, incidentally, that our approach does not assume any ad hoc postulate, contrarily to what 
believed by the author of ref . [14] . 
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At this point, wc can eventually define the mean value of the time at which our "parti- 
cle" passes through position x, travelling in the positive or negative direction of the x axis, 
respectively, as: 

< t+ix) > ^ 7^ /; \ 14.a 
/~ tJ-{x,t)dt , , 

and, moreover, the variances of the distributions of these times as: 

in accordance with the proposal presented in refs.[l,15]. 

Thus, wc have a formalism for defining mean values, variances (and other central moments) 
related to the duration distributions of all possible processes for a particle, tunnelling through a 
potential barrier located in the interval (0, a) along the x axis; and not only for tunnelling, but 
also for all possible kinds of collisions, with arbitrary energies and potentials. For instance, we 
have that 

< rT(a;i,a;f) > = < t+(xf) > - < > (16) 

with — oo < Xi < and a < Xf < oo; and therefore (as anticipated in eq.(3)) that 

DrT(xi,Xf) = Dt+(xf) + Dt+(xi) , 

for transmissions from region (— cx),0) to region (a, oo) which wc called!!] regions I and III, 
respectively. Analogously, for the pure (complete) tunnelling process one has: 



and 



while one has 



< rTun(0, a) > = < t+{a) > - < t+{0) > (17) 
DrTun(0,a) = Dt+(a) + Di+(0) ; (18) 

< rpen(0,Xf) > = < t+{xf) >-< t+{0) > (19) 
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and 



DTPe„(0,Xf) = Dt+(xf)+Dt+(0) 



(20) 



[with < Xf < a] for penetration inside the barrier region (which we called region II). Moreover: 

<TKetix,x) > = <t-{x) > - <t+{x) > (21) 



DT-R^tix,x) = Dt-{x)+Dt+{x) 



(22) 



[with < X < a] for ''''return processes''' inside the barrier. At last, for reflections in region I, 
we have that: 



< rR(xi,Xi) > = < t-{xi) > - < t+{xi) > 



(23) 



[with — oo < Xi < a], and (as anticipated in eq. (4)) that DrR(,Ti,,Ti) = Dt-{xi) + Dt+{xi). 

Let us stress that our definitions hold within the framework of conventional quantum me- 
chanics, without the introduction of any new postulates, and with the single measure expressed 
by weights (13. a), (13.6) for all time averages (both in the initial and in the final conditions). 

According to our definition, the tunnelling phase time (or, rather, the transmission du- 
ration), defined by the stationary phase approximation for quasi-monochromatic particles, is 
meaningful only in the limit x^ ^ oo when J+(x, t) is the flux density of the initial packet Jin of 
incoming waves {in absence of reflected waves, therefore). 

Analogously, the dwell time, which can be represented (cf. eqs.(l),(2)) by the 



expression 



[7,8,16] 



{Xi, X{) 



t J{xf, t) dt- t J{xi, t) dt / Jin(a;i, t) dt 

-oo J —oo J L^— oo 



with — oo < x\ < 0, and Xf > a, is not acceptable, generally speaking. In fact, the weight in the 
time averages is meaningful, positive definite and normalized to 1 only in the rare cases when 
Xi — ^ — oo and Jin = Jm (i.e., when the barrier is transparent). 



3. ~ Penetration and Return process durations, inside a rectangular barrier, for 
tunnelling gaussian wave packets: Numerical results 

We put forth here the results of our calculations of mean durations for various penetration 
(and return) processes, inside a rectangular barrier, for tunnelling gaussian wave packets; one of 
our aims being to investigate the tunnelling speeds. In our calculations, the initial wave packet 
is 

roo 

^i>i^{x,t) = G{k- k) exp[i/cx - iEt/K\ dk (24) 

^0 
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with 



G{k -k) = Cexp[-{k - kf/{2 Mf] , (25) 

exactly as in ref.[8]; and with E = fi^k'^ /2m\ quantity C being the normahzation constant, 
and m the electron mass. Our procedure of integration was described in ref.[16]. 

Let us express the penetration depth in angstroms, and the penetration time in seconds. 
In Fig. 1 we show the plots corresponding to a = 5 A, for Afe = 0.02 and 0.01 A , respectively. 
The penetration time < rpen > always tends to a saturation value. 

In Fig.2 we show, for the case Afe = 0.01 A , the plot corresponding to a = 10 A. It is 
interesting that < rpcn > is practically the same (for the same AA;) for a = 5 and a = 10 A, 
a result that confirm, let us repeat, the existence^^l of the so-called Hartman effect S"^^ Let us 
add that, when varying the parameter Ak between 0.005 and 0.15 A and letting a to assume 
values even larger than 10 A, analogous results have been always gotten. Similar calculations 
have been performed (with quite reasonable results) also for various energies E in the range 1 
to 10 eV. #2 

In Figs. 3, 4 and 5 we show the behaviour of the mean penetration and return durations 
as function of the penetration depth (with Xi = and < Xf = x < a), for barriers with 

height Vq = 10 cV and width a = 5 A or 10 A. In Fig. 3 we present the plots of < Tpcn(0, x) > 
corresponding to different values of the mean kinetic energy: E = 2.5 eV, 5 cV and 7.5 eV (plots 
1, 2 and 3, respectively) with Ak = 0.02A; and E = 5 eV with AA: = 0.04A" (plot 4), always 
with a = 5A. In Fig.4 we show the plots of < rpen(0,a;) >, corresponding to a = 5A, with 
AA; = 0.02A"^ and 0.04 A"""^ (plots 1 and 2, respectively); and to a = lOA, with AA; = 0.02A~^ 
and 0.04 A ^ (plots 3 and 4, respectively), the mean kinetic energy E being 5 eV, i.e., one half 
of Vq. In Fig. 5 the plots are shown of < T^ct{x,x) >. The curves 1, 2 and 3 correspond to 
E = 2.5 eV, 5 eV and 7.5 eV, respectively, for AA: = 0.02A and a = 5A; the curves 4, 5 
and 6 correspond to E = 2.5 eV, 5 eV and 7.5 eV, respectively, for AA; = 0.04A and a = 5A; 
while the curves 7, 8 and 9 correspond to AA: = 0.02A ^ and 0.04 A ^, respectively, ioi E = 5 
eV and a = lOA. 

Also from the new Figs. 3-5 one can see that: 1) at variance with ref.[8], no plot considered 
by us for the mean penetration duration < rpen(0, x) > of our wave packets presents any interval 
with negative values, nor with a decreasing < rpen(0, x) > for increasing x; and, moreover, 

o#2 Pqj. ^j^g interested reader, let us recall that, when integrating over dt, we used the interval 

— 10"'^'^ s to +10^^'' s (symmetrical with respect to t = 0), very much larger than the tempo- 
ral wave packet extension. [Recall that the extension in time of a wave packets is of the order of 

l/{vAk) = {Ak ^j2E/m)-^ ~ IQ-^'^ s]. Our "centroid" has been always to = 0; = 0. For 
clarity's sake, let us underline again that in our approach the initial wave packet ^'in(a:, t) is not regarded 
as prepared at a certain instant of time, but it is expected to flow through any (initial) point xi during 
the infinite time interval (— oo, -|-c»), even if with a finite time-centroid t^. The value of such centroid 
to is essentially defined by the phase of the weight amplitude G{k — k), and in our case is equal to when 
G{k — k) is real. 
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that 2) the mean tunneUing duration < rTun(0, a) > does not depend on the barrier width a 
("Hartman effect"); and finaily that 3) quantity < rTuri(0,a) > decreases when the energy 
increases. Furthermore, it is noticeabie that aiso from Figs. 3-5 we observe: 4) a rapid increase 
for the value of the electron penetration time in the initial part of the barrier region (near x = 0); 
and 5) a tendency of < rpen(0,x) > to a saturation value in the final part of the barrier, near 
X = a. 

Feature 2), firstly observed for quasi- monochromatic particles, does evidently agree with 
the predictions made in ref.[l] for arbitrary wave packets. Feature 3) is also in agreement with 
previous evaluations performed for quasi-monochromatic particles and presented, for instance, in 
refs. [1,2,15]. Features 4) and 5) can be apparently explained by interference between those initial 
penetrating and returning waves inside the barrier, whose superposition yields the resulting 
fluxes J+ and J_. In particular, if in the initial part of the barrier the returning-wave packet 
is comparatively large, it does essentially extinguish the leading edge of the incoming-wave 
packet. By contrast, if for growing x the returning-wave packet quickly vanishes, then the 
contribution of the leading edge of the incoming-wave packet to the mean penetration duration 

< ''"Pen(0, x) > does initially (quickly) grow, while in the final barrier region its increase does 
rapidly slow down. 

Furthermore, the larger is the barrier width a, the larger is the part of the back edge of the 
incoming-wave packet which is extinguished by interference with the returning-wave packet. 
Quantitatively, these phenomena will be studied elsewhere. Finally, in connection with the 
plots of < Tjiet{x,x) > as a function of x, presented in Fig.5, let us observe that: (i) the mean 
reflection duration < rR(0, 0) > = < TRct(0,0) > does not depend on the barrier width a; 
(ii) in correspondence with the barrier region betwen and approximately 0.6 a, the value of 

< TRetlO, x) > is almost constant; while (iii) its value increases with x only in the barrier region 
near x = a (even if it should be pointed out that our calculations near x = a are not so good, 
due to the very small values assumed by J_(x,t)dt therein). Let us notice that point 
(i), also observed firstly for quasi-monochromatic particles,!^] is as well in accordance with the 
results obtained in ref.[l] for arbitrary wave packets. Moreover, also points (ii) and (iii) can be 
explained by interference phenomena inside the barrier: if, near x = a, the initial returning-wave 
packet is almost totally quenched by the initial incoming— wave packet, then only a negligibly 
small piece of its back edge (consisting of the components with the smallest velocities) does 
remain. With decreasing x {x 0), the unquenched part of the returning-wave packet seems 
to become more and more large (containing more and more rapid components), thus making 
the difference < TRct(0, x) > — < rpcn(0, x) > almost constant. And the interference between 
incoming and reflected waves at points x < does effectively constitute a retarding phenomenon 
[so that < i_(x = 0) > is larger than < rR(x = 0)], which can explain the larger values of 

< T^_{x = 0, X = 0) > in comparison with < rTun(x = 0, x = a) >. 

Therefore our evaluations, in all the cases considered above, appear to confirm our previous 
analysis at page 352 of ref.[l], and our conclusions therein concerning in particular the validity 
of the Hartman effect also for non-quasi-monochromatic wave packets. Even more, since the 
interference between incoming and reflected waves before the barrier (or between penetrating 
and returning waves, inside the barrier, near the entrance wall) does just increase the tunnelling 
time as well as the transmission times, we can expect that our non-relativistic formulae for 
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< TTun(0, a) > and < TTixj < 0,Xf > a) > will always forward positive values. 



At this point, it is necessary — however — to observe the following. Even if our non- 
relativistic equations are not expected (as we have just seen) to yield negative times, neverthe- 
less one ought to bear in mind that (whenever it is met an object, O, travelling at Supcrluminal 
speed) negative contributions should be expected to the tunnelling times: and this ought not to 
be regarded as unphysical. In fact, whenever an "object" O overcomes the infinite speedt^*^! with 
respect to a certain observer, it will afterwards appear to the same observer as its "anii-object" 
O travelling in the opposite space direction!^®! . For instance, when passing from the lab to a 
frame moving in the same direction as the particles or waves entering the barrier region, the 
objects O penetrating through the final part of the barrier (with almost infinite speeds, like in 
Figs. 1-5) will appear in the frame as anti-objects O crossing that portion of the barrier in 
the opposite space-directio'd^^^ . In the new frame J^, therefore, such anti-objects O would yield 
a negative contribution to the tunnelling time: which could even result, in total, to be negative. 
For any clarifications, see refs.[18]. So, we have no objections a priori against the fact that 
Leavens can find, in certain cases, negative values^^'^^i: e.g, when applying our formulae to wave 
packets with suitable initial conditions. What we want to stress here is that the appearance of 
negative times (it being predicted by Relativity itself, ^^^^ when in presence of anything travelling 
faster than c) is not a valid reason to rule out a theoretical approach. 

At last, let us — incidentally — recall and mention the following fact. Some preliminary cal- 
culations of penetration times (inside a rectangular barrier) for tunnelling gaussian wave packets 
had been presented by us in 1994 in ref.[16]. Later on — looking for any possible explanations 
for the disagreement between the results in ref.[8] and in our ref.[16] — we discovered, however, 
that an exponential factor was missing in a term of one of the fundamental formulae on which 
the numerical computations (performed by our group in Kiev) were based: a mistake that could 
not be detected, of course, by our careful checks about the computing process. Therefore, the 
new results of ours appearing in Figs. 1-2 should replace Figs. 1-3 of ref.[16]. One may observe 
that, by using the same parameters as (or parameters very near to) the ones adopted by Leavens 
for his Figs. 3 and 4 in rcf.[8], our new, corrected figures 1 and 2 result to be more similar to 
Leavens' than the uncorrected ones (and this is of course a welcome step towards the solution 
of the problem). One can verify once more, however, that our theory appears to yield for those 
parameters non-negative results for < rpen(a;f) >, contrarily to a claim in ref.[8]. Actually, our 
previous general conclusions have not been apparently aff'ected by the mentioned mistake. In 
particular, the value of < rpen(a;f) > increases with increasing xj, and tends to saturation for 
Xf — > a. We acknowledge, however, that the difference in the adopted integration ranges [— oo 

A different claim by Delgado, Brouard and Mugal^''^ does not seem to be relevant to our calcula- 
tions, since it is based once more, like ref.[8], not on our but on different wave packets (and over-barier 
components are also retained in ref.[17], at variance with us). Moreover, in their classical example, they 
overlook the fact that the mean entrance time < i+ (0) > gets contribution mainly by the rapid compo- 
nents of the wave packet; they forget, in fact, that the slow components are (almost) totally reflected by 
the initial wall, causing a quantum-mechanical reshaping that contributes to the initial "time decrease" 
discussed by us already in the last few paragraphs of page 352 in ref.[l]. All such phenomena reduce the 
value of < t+{0) >, and we expect it to be (in our non-relativistic treatment) less than < t+{a) >. 
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to +OC for us, and to +00 for Leavens] does not play an important role, contrarily to our 
previous belief, t"*^^! in explaining the remaining discrepancy between our results and Leavens'. 
Such a discrepancy might perhaps depend on the fact that the functions to be integrated do 
fluctuate heavily*^ (anyway, we did carefully check that our own elementary integration step 
in the integration over dk was small enough in order to guarantee the stability of the numerical 
result, and, in particular, of their sign, for strongly oscillating functions in the integrand). More 
probably, the persisting disagreement can be merely due to the fact — as recently claimed also 
by Delgado et al.1^^1 — that different initial conditions for the wave packets were actually chosen 
in ref.[8] and in ref.[l]. Anyway, our approach seems to get support, at least in some particular 
cases, also by a recent article by Brouard et al., which "generalizes" — even if starting from a 
totally different point of view — some of our results. ^^^'^^1. 

Let us take advantage of the present opportunity for answering other criticisms appeared 
in ref.[8], where it has been furthermore commented about our way of performing actual av- 
erages over the physical time. We cannot agree with those comments: let us re-emphasize in 
fact that, within conventional quantum mechanics, the time t{x) at which our particle (wave 
packet) passes through the position x is "statistically distributed" with the probability densities 
dtJ±{x,t)/ J^^dtJ±{x,t), as we explained at page 350 of ref.[l]. This distribution meets the 
requirements of the time-energy uncertainty relation. 

We also answered in Sect.l Leavens' comments about our analysis^^l of the dwell-time 
approaches. [-^^1 

The last object of the criticism in ref.[8] refers to the impossibility, in our approach, of 
distinguishing between "to be transmitted" and "to be reflected" wave packets at the leading 
edge of the barrier. Actually, we do distinguish between them; only, we cannot — of course — 
separate them, due to the obvious superposition (and interference) of both wave functions in 
p{x,t), in J{x,t) and even in J±{x,t). This is known to be an unavoidable consequence of the 
superposition principle, valid for wave functions in conventional quantum mechanics. That last 
objection, therefore, should be addressed to quantum mechanics, rather then to us. Neverthe- 
less, Leavens' aim of comparing the definitions proposed by us for the tunnelling times not only 
with conventional, but also with non-standard quantum mechanics might be regarded a priori 
as stimulating and possibly worth of further investigation. 

4. — Further remarks 

In connection with the question of "causality" for relativistic tunnelling particles, let 
us stress that the Hartman-Fletcher phenomenon (very small tunnelling durations), with the 
consequence of Superluminal velocities for sufficiently wide barriers, was found theoretically 
also in QFT for Klein-Gordon and Dirac equations, and experimentally for electromagnetic 
evanescent-mode wave packetsl"^^^! (tunnelling photons). It should be recalled that the problem 
of Superluminal velocities for electromagnetic wave packets in media with anomalous dispersion, 

o#4 -y^g can only say that we succeeded in reproducing results of the type put forth in ref.[8] by using 
larger steps; whilst the "non-causal" results disappeared — in the considered cases — when adopting small 
enough integration steps. 
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with absorption, or behaving as a barrier for photons (such as regions with frustrated internal 
reflection) has been present in the scientific literature since long (sec, for instance, quotations 
[2,1,18], and refs. therein); even if a complete settlement of the causal problem is not yet avail- 
able for relativistic waves (differently from the case of point particles'^^]). Apparently, it is 
not sufficient to pay attention only to group velocity and mean duration for a particle passing 
through a medium; on the contrary, it is important taking into account and studying ab initio 
the variances (and the higher order central moments) of the duration distributions, as well as 
the wave packet reshaping in presence of a barrier, or inside anomalous media (even if reshaping 
does not play always an essential role) . 

Passing to the approaches alternative to the direct description of tunnelling processes in 
terms of wave packets, let us here recall those ones which are based on averaging over the set of 
all dynamical paths (through the Feynman path integral formulation, the Wigner distribution 
method, and the non-conventional Bohm approach), and others that use additional degrees of 
freedom which can be used as "clocks" . General analyses of all such alternative approaches can 
be found in refs.[l, 20-24] from different points of view. 

If one confines himself within the framework of conventional quantum mechanics, then the 
Feynman path integral formulation seems to be adequate, l^^l But it is not clear what procedure 
is needed to calculate physical quantities within the Feynman-type approach l^^l, and usually 
such calculations result in complex tunnelling durations. The Feynman approach seems to need 
further modifications if one wants to apply it to the time analysis of tunnelling processes, and 
its results obtaind up to now cannot be considered as final. 

As to the approaches based on introducing additional degrees of freedom as "clocks", one 
can often realize that the tunnelling time happens to be noticeably distorted by the presence of 
such degrees of freedom. For example, the Buttiker-Landauer time is connected with absorption 
or emission of modulation quanta (caused by the time-dependent oscillating part of the barrier 
potential) during tunnelling, rather than with the tunnelling process itself. And, with 

reference to the Larmor precession time, it has been shownl^^'^'^l that this time definition is 
connected not only with the intrinsic tunnelling process, but also with the geometric boundaries 
of the magnetic field introduced as a part of the clock: for instance, if the magnetic field region 
is infinite, one ends up with the phase tunnelling time, after an average over the (small) energy 
spread of the wave packet. Actually, those "clock" approaches, when applied to tunnelling 
wave packets, seem to lead — after eliminating the distortion caused by the additional degrees 
of freedom — to the same results as the direct wave packet approach, whatever be the weight 
function adopted in the time integration. 
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Figure Captions 



Fig.l - Behaviour of the average penetration time < Tpen(0, x) > (expressed in seconds) as a 
function of the penetration depth Xi = x (expressed in angstroms) through a rectangular barrier 
with width a = 5 A, for Ak = 0.02 A (dashed Unc) and Ak = 0.01 A (contimious Hne), 
respectively. The other parameters are listed in footnote ^1. It is worthwhile to notice that 
< Tpen > rapidly increases for the first, few initial angstroms (~ 2.5 A), tending afterwards to 
a saturation value. This seems to confirm the existence of the so-called "Hartman effect". I^'^'^^l 

Fig.2 - The same plot as in Fig.l, for AA; = 0.01 A ^, except that now the barrier width 
is a = 10 A. Let us observe that the numerical values of the (total) tunnelling time < tt > 
practically does not change when passing from a = 5 A to a = 10 A, again in agreement with 
the characteristic features^ of the Hartman effect. Figures 1 and 2 do improve (and correct) 
the corresponding ones, preliminarily presented by us in ref.[16]. 

Fig.3 - Behaviour of < Tpen(0,x) > (expressed in seconds) as a function of x (expressed in 
angstroms), for tunnelling through a rectangular barrier with width a = 5 A and for different 
values of E and of Afc: 

curve 1: Ak = 0.02A"^ and E = 2.5 eV; curve 2: AA; = 0.02A"^ and E = 5.0 eV; curve 3: 
Ak = 0.02A"^ and E = 7.5 eV; curve 4: Ak = 0.04A"^ and E = 5.0 eV. 

Fig.4 - Behaviour of < rpen(0, x) > (in seconds) as a function of x (in angstroms) for £^ = 5 eV 
and different values of a and Ak: 

curve 1: a = 5 A and AA; = 0.02A~^; curve 2: a = 5 A and AA; = 0.04A~'^; curve 3: a = 10 A 
and AA = 0.02A"^; curve 4: a = 10 A and AA = 0.04A"\ 

Fig. 5 -Behaviour of < T^et{x,x) > (in seconds) as a function of x (in angstroms) for different 
values of a, E and AA;: 

curve 1: a = 5 A, E = 2.5 eV and AA = 0.02A~^ curve 2: a = 5 A, E = 5.0 cV and 
AA = 0.02A^^; curve 3: a = 5 A, E = 7.5 eV and AA = 0.02A^^; curve 4: a = 5 A, 
E = 2.5 eV and AA = 0.02A"^ curve 5: a = 5 A, E = 5.0 eV and AA = 0.02A"^ curve 6: 
a = 5 A, E = 7.5eV and AA = 0.02A"^ curve 7: a = 5 A, E = 5.0 eV and AA = 0.02A"^ 
curve 8: a = 5 A, E = 5.0 eV and AA; = 0.02A~\ 
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